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1 Introduction 

Turbulence is one of the most interesting research fields in mechanics. But with the current per- 
formance of computers, a direct simulation of a turbulent flow remains difficult, even impossible 
in many cases, due to the high computational cost that it requires. To reduce this computational 
cost, use of turbulence models is necessary. At the present time, many turbulence models exist 
(see |22j). However, derivation of a very large majority of them does not take into account the 
symmetry group of the basic equations, the Navier-Stokes equations. 

In turbulence, symmetries play a fundamental role in the description of the physics of the 
flow. They reflect existence of conservation laws, via Noether's theorem. Notice that, even if 
Navier-Stokes equations are not directly derived from a Lagrangian, Noether's theorem can be 
applied and conservation laws can be deduced. Indeed, there exists a Lagrangian (which will be 
called "bi-Lagrangian" here) from which Navier-Stokes equations, associated to their "adjoint" 
equations, can be derived. A way in which this bi-Lagrangian can be calculated is described 
by Atherton and Homsy in pQ. The expression of this bi-Lagrangian and the infinitesimal 
generators of the associated Euler-Lagrange equations are given in Appendix A. However, the 
conservation laws are not studied in this paper. 

The importance of symmetries in turbulence is not limited to the derivation of conservation 
laws. Unal also used a symmetry approach to show that Navier-Stokes equations may have 
solutions which have the Kolmogorov form of the energy spectrum [21]. Next, symmetries 
enabled Oberlack to derive some scaling laws for the velocity and the two point correlations |19j . 
Some of these scaling laws was reused by Lindgren et al in |15j and are proved to be in good 
agreement with experimental data. Next, symmetries allowed Fushchych and Popowych to 
obtain analytical solutions of Navier-Stokes equations [S]. The study of self-similar solutions 
gives also an information on the behaviour of the flow at a large time (I]. Lastly, we mention 
that use of discretisation schemes which are compatible with the symmetries of an equation 
reduces the numerical errors |2(J1 113) . 
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Introduction of a turbulence model in Navier-Stokes equations may destroy symmetry pro- 
perties of the equations. In this case, physical properties (conservation laws, scaling laws, 
spectral properties, large-time behaviour, . . . ) may be lost. In order to represent the flow 
correctly, turbulence models should then preserve the symmetries of Navier-Stokes equations. 
The first aim of this paper is to show that most of the commonly used subgrid turbulence 
models do not have this property. The second goal is to present a new way of deriving models 
which are compatible with the symmetries of Navier-Stokes equations and which, unlike many 
existing models, conform to the second law of thermodynamics. As it will be shown in appendix, 
conformity with this law leads to stability of the model, in the sense of L 2 . 

The paper will be structured as follows. In Section 2, the principle of turbulence modelling, 
using the large-eddy simulation approach, will be concisely presented, as well as some common 
models. These models will be analysed in Section 3 under the symmetry consideration. In 
Section 4, a class of symmetry preserving and thermodynamically consistent models are derived. 
One example of a model of the class is numerically tested in Section 5. Some conclusions will 
be drawn in Section 6. In Appendix A, it will be shown that Navier-Stokes equations can be 
derived from a bi-Lagrangian. At last, in Appendix B, stability of thermodynamically consistent 
models is proved. 



2 Large-eddy simulation 

Consider a three-dimensional incompressible Newtonian fluid, with density p and kinematic 
viscosity v. The motion of this fluid is governed by Navier-Stokes equations: 

du 1 

— — h divftt (g) u) H — Vp = div T, 

ot p 

divu = 0, (1) 

where u = {ui)i=i^,z and p are respectively velocity and pressure fields and t the time variable. 
T is a tensor such that pT is the viscous constraint tensor. T can be linked to the strain rate 
tensor S = (Vu + T Vit) /2 according to the relation: 

tp being a positive and convex "potential" defined by: 
ip = vti S 2 . 



Since a direct numerical simulation of a realistic fluid flow requires very significant compu- 
tational cost, (^Q) is not directly resolved. To circumvent the problem, some methods exist. The 
most promising one is the large-eddy simulation. It consists in representing only the large scales 
of the flow. Small scales are dropped from the simulation; however, their effects on the large 
scales are taken into account. This enables to take a much coarser grid. 

Mathematically, dropping small scales means applying a low-pass filter. The large or resolved 
scales of a quantity <j) are defined by the convolution: 

4> = Gg * <j), 

where Gj is the filter kernel with a width 5, and the small scales 0' are defined by 



<// = </>-<!>. 
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It is required that the integral of Gj over ]R 3 is equal to 1, such that a constant remains unchanged 
when the filter is applied. 

In practice, (u,p) is directly used as an approximation of (u,p). To obtain (u,p), the filter 
is applied to If the filter is assumed to commute with the derivative operators (that is not 
always the case in a bounded domain), this leads to: 

du 1 _ 

— + div(u ® u) + -Vp = div(T + T s ), 

at p 

divu = 0, (2) 



where T s is the subgrid stress tensor defined by T s = u ®u — u® u which must be model- 
led (expressed by a function of the resolved quantities) to close the equations. Currently, an 
important number of models exists. Some of the most common ones will be reminded here. 
They will be classified in four categories: turbulent viscosity, gradient-type, similarity- type and 
Lund-Novikov-type models. 

2.1 Turbulent viscosity models 

Turbulent viscosity models are models which can be written in the following form: 

Tf = u s S, (3) 
where v s is the turbulent viscosity. The superscript ( d ) represents the deviatoric part of a tensor: 

Q ^ Q rf = Q-I(trQ)l d , 

where Id is the identity operator. The deviatoric part has been introduced in order to have the 
equality of the traces in ©. In what follows, some examples of turbulent viscosity models are 
presented. 

• Smagorinsky model (see |22j ) is one of the most widely used models. It uses the local 
equilibrium hypothesis for the calculation of the turbulent viscosity. It has the following 
expression: 

Tf = (C S 5) 2 \S\S, 



where Cs — 0.148 is the Smagorinsky constant, 5 the filter width and |S| = y 2tr (S 2 ). 

In order to reduce the modelling error of Smagorinsky model, Lilly ^1] proposes a dynamic 
evaluation of the constant Cs by a least-square approach. This leads to the so-called 
dynamic model defined by: 

a -2-- tr(LM) 
Tf = CRISIS, with Q = -A_Z (4) 

In these terms, 

L = u®u-u®«, M = (5 2 1 S | S — 5 |S|S, 

the tilde represents a test filter whose width is 5, with 5 > 5. 

The last turbulent viscosity model which will be considered is the structure function model. 
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• Metais and Lesieur JS] make the hypothesis that the turbulent viscosity depends on the 
energy at the cutoff. Knowing its relation with the energy density in Fourier space, they 
use the second order structure function and propose the structure function model: 

T d s = C SF 5^J¥^)S, (5) 

where F2 is the spatial average of the filtered structure function: 

r 1 — ^ F<z(r) = // \\u(x) — u(x + z)\\ 2 dz dx. 

J J \\z\\=r 



The next category of models, which will be reminded, consists of the gradient-type models. 



2.2 Gradient-type models 

To establish the gradient-type models, the subgrid stress tensor is decomposed as follows: 

T s = tt ® w — (u®u + u®u' + u'®u + u'<giu'). 

Next, each term between the brackets are written in Fourier space. Then, the Fourier transform 
of the filter, which is assumed to be Gaussian, is approximated by an appropriate function. 
Finally, the inverse Fourier transform is computed. The models in this category differ by the 
way in which the Fourier transform of the filter is approximated. 

• If a second order Taylor series expansions according to the filter width 5 is used in the 
approximation, one has: 

5 2 

T s = -—Vu T Vu. (6) 



• The gradient model is not dissipative enough and not numerically stable |26Mllj . Thus, it 
is generally combined to Smagorinsky model. This gives Taylor model: 

■f 

T s = Vu T VTZ+C5 2 |S|S. 



The Taylor approximation of the Fourier transform of the filter tends to accentuate the 
small frequencies rather than attenuating them. Instead, a rational approximation can be 
used [2] . This gives the following expression of the model: 

5 2 ( 5 2 \ l 
T s = --I l d --V 2 j [Vu T Vu]+cf\S\S. 

To avoid the inversion of the operator ^Ij — I^V 2 ^, T s is approximated by: 
5 2 

T s = -—G I *[Vu T Vu]+C5 2 \S\S. (7) 

G-g is the kernel of the Gaussian filter. The convolution is done numerically. The model 
is called the rational model. 
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2.3 Similarity-type models 

Models of this category are based on the hypothesis that the statistic structure of the small 
scales are similar to the statistic structure of the smallest resolved scales. Separation of the 
resolved scales is done using a test filter (symbolized by ~). The largest resolved scales are then 
represented by u and the smallest ones by u — u. From this hypothesis, we deduce the similarity 
model: 

T s = u®u — u®u. (8) 

From this expression, many other models can be obtained by multiplying by a coefficient, by 
filtering again the whole expression or by mixing with a Smagorinsky-type model. 
The last models that we will consider are Lund-Novikov-type models. 

2.4 Lund— Novikov-type models 

• Lund and Novikov include the filtered vorticity tensor W = (Vm — T ViZ) in the expres- 
sion of the subgrid model. Cay ley-Hamilton theorem gives then the Lund-Novikov model 



where the coefficients C{ depend on the invariants obtained from S and W. The expression 
of these coefficients are so complex that they are considered as constants and evaluated 
with statistic techniques. 

• To reduce the computation cost of the previous model, Kosovic brings a simplification and 
proposes the following model: 



where the constants C, C\ and C2 are calculated using the theory of homogeneous and 
isotropic turbulence. 

The derivation of these models was done using different hypothesis but did not take into 
consideration the symmetries of Navier-Stokes equations which may then be destroyed. So, in 
the next section, these models will be analysed by a symmetry approach. 

3 Model analysis 

The (classical) symmetry groups of Navier-Stokes equations have been investigated for some 
decades (see for example QUI]). They are generated by the following transformations: 

• The time translations: (t,x,u,p) \—t (t + a,x,u,p), 

• the pressure translations: (t,x,u,p) 1— > (t,x,u,p + £(t)), 

• the rotations: (t, x, u,p) 1— > (t, Rx, Ru,p), 

• the generalized Galilean transformations: 

(t, x,u,p) i — > (t, x + oc(t), u + ac(t), p — px ■ ac(t)), 

• and the first scaling transformations: (t,x,u,p) 1— > (e 2a t, e a x, e~ a u, e~ 2a p). 



(see |22|): 



T d s = d5 2 \S\S + C 2 <5 2 (S 2 ) d + C 3 5 2 (\N 2 ) d 




(9) 



T d s = (CS) 2 [2|S|S + Ci(S 2 ) d + C 2 (SW- WS) 



(10) 
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In these expressions, a is a scalar, £ (respectively a) a scalar (resp. vectorial) arbitrary 
function of t and R a rotation matrix, i.e. R T R = and det R = 1. The central dot (•) stands 
for M 3 scalar product. 

If it is considered that v can change during the transformation (which is then an equivalence 
transformation [5]), one has the second scaling transformations: 

(t,x,u,p,v) i ^ (t,e a x,e a u,e 2a p,e 2a u), 

where a is the parameter. 

Navier-Stokes equations admit other known symmetries which do not constitute a one- 
parameter symmetry group. They are 

• the reflections: (t, x, u,p) i— > (t, /\x, Ait,p), 

which are discrete symmetries, A being a diagonal matrix A = diag(ti, t2 ; ^3) with Li = 
±1, i = 1,2,3, 

• and the material indifference: (t,x,u,p) 1— > (t,x,u,p), 
in the limit of a 2D flow in a simply connected domain with 

x = R(t) x, u = R(t) u + R(t) x, p = p — 3uj(p + -^u 2 \\x\\ 2 , 

where R(t) is a 2D rotation matrix with angle tot, to an arbitrary real constant, <p the usual 
2D stream function defined by: 

u = curl((/?e3), 

63 the unit vector perpendicular to the plane of the flow and ||«|| the Euclidean norm. 

We wish to analyse which of the models cited above is compatible with these symmetries. 
The set of solutions (u,p) of Navier-Stokes equations (^Q) is preserved by each of the symmetries. 
We then require that the set of solutions (u, p) of the filtered equations ((2j) is also preserved by 
all of these transformations, since (u,p) is expected to be a good approximation of (u,p). More 
clearly, if a transformation 

T : (t,x,u,p) 1— > (t,x,u,p) 

is a symmetry of we require that the model is such that the same transformation, applied 
to the filtered quantities: 

T : (t,x,u,p) 1 ^ (t,x,u,p), 

is a symmetry of the filtered equations @- When this condition holds, the model will be said 
invariant under the relevant symmetry. 

The filtered equations © may have other symmetries but with the above requirement, we 
may expect to preserve certain properties of Navier-Stokes equations (conservation laws, wall 
laws, exact solutions, spectra properties, . . . ) when approximating (u,p) by (It, p). 

We will use the hypothesis that test filters do not destroy symmetry properties, i.e. 4> = 4> 
for any quantity (j>. 

For the analysis, the symmetries of will be grouped into four categories: 

- translations, containing time translations, pressure translations and the generalized Gali- 
lean transformations, 
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- rotations and reflections, 

- scaling transformations, 

- material indifference. 

The aim is to search which models are invariant under the symmetries within the considered 
category. 

3.1 Invariance under translations 

Since almost all existing models are autonomous in time and pressure, the filtered equations (J2J) 
remain unchanged when a time or pressure translation is applied. Almost all models are then 
invariant under the time and the pressure translations. 

The generalized Galilean transformations, applied to the filtered variables, have the following 



All models in Section 2, in which x and u are present only through Vu. are invariant since 



where V = (8/8x1,8/8x2,6/8x3). 

The remaining models, i.e. the dynamic and the similarity models are also invariant because 

u(g>tt — u ®u = (u + a) ® (u + a) — (u + a) ® (u + a) = u ® u — u ® u. 

3.2 Invariance under rotations and reflections 

The rotations and the reflections can be put together in a transformation: 

(t, x, u,p) 1 — ^ (t, Tx, Tu,p) 

where T is a constant rotation or reflection matrix. This transformation, when applied to the 
filtered variables, is a symmetry of J2J) if and only if 



form: 



(t,x,u,p) 1 ► (t,x,u,p) = (t,x + a(t),u + at(t),p — px • ix(t)) . 



Vu = Vu = Vu, 



T T/T. 



(11) 



Let us check if the models respect this condition. 



• For Smagorinsky model, we have: 



Vu = [V(u)] 



T T = [V(Tw)] T T = T[ViZ] T T. 



(12) 



This leads to the objectivity of S: 



S = T S T T. 



And since IS 



S|, Ijllj) is verified. Smagorinsky model is then invariant. 



• For similarity model ©, one has: 



u <g> u = (Tit) (g> (Tm) = T(u £$ u) 



By means of these relations, invariance can easily been deduced. 
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• The same relations are sufficient to prove invariance of the dynamic model since the trace 
remains invariant under a change of orthonormal basis. 

• The structure function model © is invariant because the function F2 is not altered under 
a rotation or a reflection. 

• Relations (|12|) can be used again to prove invariance of each of the gradient-type models. 

• Finally, since 

W= T W ' T. 

Lund-Novikov-type models are also invariant. 
Any model of Section 2 is then invariant under the rotations and the reflections. 

3.3 Invariance under scaling transformations 

The two scaling transformations can be gathered in a two-parameter transformation which, when 
applied to the filtered variables, have the following expression: 

{t,x,u,p,u) 1 ^ (e 2 %e ab x,e b - a u,e 2b - 2a p,e 2b v). 

where a and b are the parameters. The first scaling transformations corresponds to the case 
6 = and the second ones to the case a = 0. 

It can be checked that the filtered equations (J2J) are invariant under the two scaling trans- 
formations if and only if 

f s = e 2b ~ 2a T s . (13) 

Since S = e _2a S, this condition is equivalent to: 

v s = e 2b u s (14) 
for a turbulent viscosity model. 

• For Smagorinsky model, we have: 

v s = C s f\%\ = e- 2a C s f\S\ = e" 2 V s . 

Condition (|14j) is violated. The model is invariant neither under the first nor under the 
second scaling transformations. Note that the filter width 5 does not vary since it is an 
external scale length and has no functional dependence on the variables of the flow. 

• The dynamic procedure used in (jlj) restores the scaling invariance. Indeed, it can be shown 
that: 

C d = e 2b+2a C d , 
that implies: 

u s = C d 5 2 \i\ = e 2b C d 5 2 \S\ = e 2b v s . 
The dynamic model is then invariant under the two scaling transformation. 
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• For the structure function model, we have: 

f 2 = e b ~ a F 2 
and then 

V s — e fsmi 

that proves that the model is not invariant. 

• Since 

VSZ = e 2a Vu 

the gradient model (jHJ) violates 1)13(1. T s varying in the following way: 
f s = e 4a T s . 

This also implies that none of the gradient-type models is invariant. 

• It is straight forward to prove that the similarity model © verifies (jl3|) and is invariant. 

• At last, Lund-Novikov-type models are not invariant because they comprise a term similar 
to Smagorinsky model. 

In fact, none of the models where the external length scale 8 appears explicitly is invariant 
under the scaling transformations. Note that the dynamic model, which is invariant under these 
transformations, can be written in the following form: 

, tr(LN) __ 
j d = : L SIS 

s tr(N2) 

where 

n = jiji- (f/5) 2 |l|f. 

It is then the ratio 8/8 which is present in the model but neither 8 alone nor 8 alone. 

In summary, the dynamic and the similarity models are the only invariant models under the 
scaling transformations. Though, scaling transformations have a particular importance because 
it is with these symmetries that Oberlack derived scaling laws and that Unal j21] proved 
the existence of solutions of Navier-Stokes equations having Kolmogorov spectrum. 

The last symmetry property of Navier-Stokes equations is the material indifference, in the 
limit of 2D flow, in a simply connected domain. 



3.4 Material indifference 

The material indifference corresponds to a time-dependent plane rotation, with a compensation 
in the pressure term. We will not write explicitly the dependence on time of the rotation 
matrix R. 

• The objectivity of S (see Section l3~2() directly leads to invariance of Smagorinsky model. 
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• For similarity model (JHJ), we have: 

f s = RT/R + R(ll®rx - M <g) x) T R 

+ R(x tg) u - 2 <g> u) T R + R(x ~®x — x® x) T R. 

Consequently, if the test filter is such that 

(u <g> x — u g> x) = 0, (x (g u — x (g w) = 0, (a? (g a; — 5 <g> a;) = 0, (15) 

then the similarity model is invariant under the material indifference. All filters do not 
have this property. For instance, it can be shown |2j that, for the usual box filter, the 

~ ~ ~2 

left-hand sides of equations (|15|) do not vanish and are respectively in 0(5), 0(5) and 0(5 ). 

• Under the same conditions ()15|) on the test filter, the dynamic model is also invariant. 

• The structure function model is invariant if and only if 

f 2 = F 2 . (16) 
Let us calculate F 2 . Let u z be the function x i— > u(x + z). Then 
%= _\\(Ru + Rx) - (Ru z + Rx + Rz)\\ 2 dz 

J\\z\\=5 

= / ||tt — u z — T RR,z|| dz. 

J\\z\\=5 

Knowing that T RRz = uje-s x z, we get: 

F 2 = F 2 + 2ituj 2 5 3 - 2u I _(u - u z ) ■ (e 3 x z) dz. 

J\\z\\=S 

Condition (|16j) is violated. So, the structure function model is not invariant under the 
material indifference. 

• For the gradient model, we have: 

W = RVU T R + R T R, T W= R T Vw T R + R T R. (17) 
Let J be the matrix such that R T R = — ooS = — R T R or, in a component form: 

V-l 

Then, 

(Vw^Vti) = R(Vm t Vu) t R + cuRVm t RJ - wJR T Vw T R + u 2 l d . 
The commutativity between J and R finally leads to: 
f s = R T s T R + wR(VwJ - J T Vu) + uj 2 \ d . 

This proves that the gradient model is not invariant. 

• The other gradient-type models inherit the lack of invariance of the gradient model. 
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• It remains the Lund-Novikov-type models. We will begin with Kosovic model (|ll)j) since 
it is simpler. The first two terms of (|1U|) are unchanged under the transformation. For the 
filtered vorticity tensor W, it follows from Q17JI that: 



W = RW T R-wJ. 



Thus, 



SW-WS = R(SW-W S) T R - wR(SJ - JS) T R 

using again the commutativity between J and R. As for them, S and J are not com- 
mutative. In fact, using properties of S, it can be shown that SJ = — JS. This implies 
that 



SW-WS = R(SW-W S) T R - 2cuRSJ T R. 

This shows that Kosovic model is not invariant. 
Lastly, consider Lund-Novikov model @. We have: 

W 2 = RW 2T R-wR(jW + Wj) T R-u; 2 l d . 

Since W is anti-symmetric and the flow is 2D, W is in the form: 

W =( - n 

\—w 

A direct calculation leads then to 
JW = WJ = -wl d 

and 



(18) 



W 



RW T R- (2w-u)u\ a . 



Let us see now how each term of the model © containing W varies under the transfor- 
mation. 

From the last equation, we deduce the objectivity of (\A/ 2 ) d : 

(w 2 )^ = r (wy t r. 

For the fourth term of Q, we already have (|18j). And for the last term, 



SW-SW = R(S'W- SW ) T R -cjRS T RJ - (2w - lu)ujRS t R. 

Putting these results together, we have: 



RT d a t R-uj5 R 



2C 4 SJ - C 5l =r ( S'J - (2w - oo)ujS 



R. 



We conclude that Lund-Novikov model is not invariant the material indifference. This 
ends the analysis. 
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Translations 


Rotations, 
reflections 


Scaling 
transformations 


Material 
indifference 


Srnagorinsky 


Y 


Y 


TV 


Y 


Dynamic 


Y 


Y 


Y 




Structure function 


Y 


Y 


TV 


N 


Gradient 


Y 


Y 


TV 


N 


Taylor 


Y 


Y 


N 


N 


Rational 


Y 


Y 


N 


N 


Similarity 


Y 


Y 


Y 


•y* 


Lund 


Y 


Y 


N 


N 


Kosovic 


Y 


Y 


N 


N 



Table 1. Results of the model analysis. 
Y=invariant. 7V=not invariant, Y*=invariant if (|15|) is verified. 

Table n summarizes the results of the above analysis. "Y" means that the model is invariant 
under all the symmetries of the category, "JV" the opposite and " Y*" that the model is invariant 
if the conditions (|15j) on the test filter is verified. It can be seen on this table that only two models 
among the nine, the dynamic and the similarity models, are invariant under the symmetry group 
of Navier-Stokes equations. The scaling transformations, which are of a particular importance 
(scaling laws, Kolmogorov spectrum, . . . ) are violated by almost all models. 

The dynamic and the similarity models have an inconvenience that they necessitate use of a 
test filter. Rather constraining conditions, 1)15(1 . are then needed for these models to preserve 
the material indifference. In addition, the dynamic model does not conform to the second law of 
thermodynamics since it may induce a negative dissipation. Indeed, v + v s can take a negative 
value. To avoid it, an a posteriori forcing is generally done. It consists of assigning to u s a value 
slightly higher than — v: 

v s = -v (1 - e), 

where e is a positive real number, small against 1. Non-conformity to the second law of thermo- 
dynamics may be detrimental for a model because, as it will be shown in Appendix B, consistence 
with this law leads to stability of the model. 

Considering this lack of invariance of existing models and to non-conformity with thermo- 
dynamical principles, we propose in the next section a new way of deriving models which, on 
one hand, possess the symmetry group of Navier-Stokes equations and, on the other hand, are 
compatible with the second law of thermodynamics. 

4 Invariant and thermodynamically consistent models 

First, we will build a class of models which possess the symmetries of Navier-Stokes equations 
and next refine this class such that the models also satisfy the thermodynamics requirement. 
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4.1 Invariance under the symmetries 

Suppose that S / 0. Let T s be an analytic function of S: 

T S =.4(S). (19) 

By this way, invariance under the time, pressure and generalised Galilean translations and under 
the reflections is guaranteed. Prom (|19|). Cayley-Hamilton theorem and invariance under the 
rotations lead to: 

Tf = A(xX)S + B(x,()Mj d S, (20) 

where \ = trS 2 and £ = detS are the invariants of S (the third invariant, trS, vanishes), Adj 
stands for the operator defined by 

(AdjS)S = (detS)ld, 

(Adj S is simply the comatrix of S ) and A and B are arbitrary scalar functions. Contrarily to 
Lund-Novikov model, these coefficient functions will not be taken constant. 

Next, a necessary and sufficient condition for T s defined by (|2U|) to be invariant under the 
second scale transformations is that v can be factorized: 

Tf = uAo(x, C) S + vB ( X , C) Adj d 5. 
Lastly, T s is invariant under the first scaling transformations if 
f s = e~ 2a T s . 
Rewritten for Aq and Bq, this condition becomes: 

A (e- 4a X , e~ 6a = Mx, 0, A)^*, e" 6a C) = e 2a B ( X , ()■ 

After differentiating according to a and taking a = 0, it follows: 

dA dB dB 

ox d( ox o( 

To satisfy these equalities, one can take 

MX,0=A 1 I-^), B (x,() = -^B 1 



x 3/2 I ' ^ I x 3/2 



Finally, if v = ^, then 

v / 



Tf = vAAv) S + v—B^A&fS. 

Vx 



(21) 



A subgrid-scale model of class ()21[) remains then invariant under the symmetry transforma- 
tions of Navier-Stokes equations. 

In fact several authors were interested in building invariant models for a long time. But 
because they did not use Lie theory, they did not consider some symmetries such as the sca- 
ling transformations which are particularly important. Three of the few authors who consi- 
dered all the above symmetries in the modeling of turbulence are Unal and Saveliev and 
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Gorokhovski [23]. The present manner to build invariant models generalises the Unal's one in 
the sense that it introduces v and the invariants of S into the models. In addition, Unal used the 
Reynolds averaging approach (RANS) instead of the large-eddy simulation approach (LES) for 
the turbulence modelling. Saveliev and Gorokhovski in |23| used the LES approach but derive 
their model in a different way than in the present article. 

Let us now return to considerations which are more specific to large eddy simulation. We 
know that T s represents the energy exchange between the resolved and the subgrid scales. 
Then, it generates certain dissipation. To account for the second law of thermodynamics, we 
must ensure that the total dissipation remains positive that is not always verified by models in 
the literature. In order to satisfy this condition, we refine class (|21|). 



4.2 Consequences of the second law of thermodynamics 

At molecular scale, the viscous constraint is: 

dib 

The potential if) = v tr S 2 is convex and positive that ensures that the molecular dissipation is 
positive: 

$ = tr(TS) > 0. 

The tensor T s can be considered as a subgrid constraint, generating a dissipation 
$ s = tr(T s S). 

To preserve compatibility with the Navier-Stokes equations, we assume that T s has the same 
form as T: 

T s = ^i. (22) 

where ip s is a potential depending on the invariants x an d C °f S. This hypothesis refines 
class (j21|) in the following way. 

Since trS = 0, one deduces from (f2*2*|) : 

Comparing it with (|21|) . one gets: 

1 dt/j s 1 dij} s 
-vA l (v) = —, v—Bi{v) = —. 

2 dx y/x d( 

This leads to: 

I (^w) -I (><»>)■ 

If g is a primitive of B\ , a solution of this equation is 



Ai(u) = 2g(v) - 3vg'(v) and B^v) = g'(v). 



(23) 
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Then, the hypothesis (|22|) involves existence of a function g such that: 

(24) 



T d s = v [2g(v) - 3vg'(v)} S + u— <?» Adj d S. 



Now, let $r be the total dissipation. We have: 
$ T = tr[(T + T s )S]. 
Using (|21j) et (|23ft . one can show that 



<S> T > <^ 1 + Ai(u) + 3vSi(u) > 

1 + > 0. (25) 

In summary, a model belonging to class (|24|) with a continuous function g verifying 



l + g>0, (26) 



is a model possessing the symmetry group of Navier-Stokes equations conform to the second law 
of thermodynamics. Such a model can take into account the inverse energy cascade, since & s 
can have negative values. In addition, by putting S in a diagonal form, it can be shown that 
v belongs to a bounded interval [— v*,v*] where v* ~ 0.136. Consequently, it is not necessary 
to satisfy (|26|) out of this interval. Another important property of such a model is its stability, 
in the sense that the L 2 -norm of the filtered velocity remains bounded. In fact, all models 
which are consistent with the second law of thermodynamics are stable. This will be proved in 
Appendix B. 

In the next section, we show that our approach can lead to numerically efficient results. 
A very simple model of class Q24|) is then chosen and compared to the two most popular models 
which are Smagorinsky and the dynamic models (see |22| I14j ) . 



5 Numerical test 

We choose a simple linear function for g: 
g (v) = Cv. 

where the constant C can depend on the filter width and other parameters. 
Let d be the ratio: 




where I is a length scale related to the size of the domain. The introduction of this ratio is also 
useful to have the right dimensions. We now take: 

C = (C s d) 2 , 

where C s is a pure constant, set to be equal to Smagorinsky constant, i.e. C s ~ 0.16. Doing so, 
condition (|25|) is verified and one has: 
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L=9m 



inlet velocity : U jn = 0.455 m/s 
inlet height : 0.168 m 
outlet height : 0.46 m 





-0.4 



-0.2 
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0.4 



0.6 



Velocity ;/j AJj t 



m 



Figure 1. Geometry of the ventilated room. 



Figure 2. Mean velocity profiles at x/L = 2/3. 



In the present paper, this model will be called "invariant model". We use this model to 
simulate a flow within a ventilated room (Nielsen's cavity JZ]) which interests us particularly 
for applications in building field. The results will then be compared to those provided by the 
Smagorinsky model and the dynamic model. 

The geometry of the room is presented on Fig. ^ For this configuration, we take £ = lm. 

The code used for the resolution was developed by Chen et al and is described in f]]. The 
spatial discretization is performed by a finite difference scheme. 

Fig. 13 compares the velocity profiles given by Smagorinsky, the dynamic and the invariant 
models with experimental data at x±/L = 2/3 and x%/W = 0.5. It can be observed on it that 
the invariant model gives a better result than Smagorinsky and dynamic models, without need 
of a test filtering. The result is in good agreement with experiments, except near the floor. 
Notice that no wall model was used. 

6 Conclusion 

In this article, we presented a new class of physically compatible subgrid turbulence models. The 
main ingredient used is the symmetry group of the Navier-Stokes equations which contains a 
fundamental information on the properties of the flow. The second principle of thermodynamics 
was also introduced. From a practical point of view, conformity with this principle ensures 
stability of the model. 

A simple model of the class was tested and encouraging results was obtained. However, the 
aim of this test was not to present a complete analysis of the model but to check that the 
symmetry approach can lead to good numerical results. Further studies will be done in future 
works on the choice of the parameters of the model and on analysis of the numerical results. 

The way presented here for deriving symmetry compatible models is a general way. It can be 
applied to other equations (non-isothermal fluid, . . . ). Other parameters can also be included. 
For example, dependence of the model on the viscosity v can be replaced by dependence on the 
dissipation rate. 

A Noether's theorem to Navier— Stokes equations 

Noether's theorem can be applied to evolution equations which can be derived from a Lagrangian, 
i.e. evolution equations which can be expressed in an Euler-Lagrange form: 




Or 



dr 



(27) 
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where r is the dependent variable, r = r(y), y = (yi)i the independent variable, L the La- 
grangian and Div the operator: 

/^Div/ = V-i. 

From the infinitesimal generators of (|27|). conservation laws are deduced. 

Navier-Stokes equations cannot be directly written in the form (|27|). However, thanks to an 
approach of Atherton and Homsy pQ , see also JO] i which consists in extending the Lagrangian 
notion, it will be shown in this appendix that Noether's theorem can be applied to Navier-Stokes 
equations. 

We will say that an evolution equation 

F(y,r)=0 (28) 
is derived from a "bi-Lagrangian" if there exists an (non necessarily unique) application 

L : (r,s) i — ^ L(r,s) G R 

such that 1)28|) is equivalent to 

9L(r,a) dL(r,s) 

a Dlv 7T- = °- 

OS OS 

s is called the adjoint variable and the equation 

dL(r,s) _ dL{r,s) = 
dr dr 

is called the adjoint equation of (|28|). The Noether theorem can then be applied since the 
evolution equation, associated to his adjoint, can be written in an Euler-Lagrangian form 

dL(w) dL(w) 

— Div . = 0, 

ow ow 

where w = (r, s). 

Navier-Stokes equations are derived from a bi-Lagrangian 

1 ( du dv\ ( 1 \ 

L((u,p), ( V >V)) = ' v ~ u ' ~dt ) + [ q ~ 2 U ' V ) ~ pd[vv + ^ tr ( TVu • Vtl ) • 

where i> = (fj), and q are the adjoint variables. The corresponding adjoint equations are 
dv 

-— + (v ■ Vu T - u ■ Vv) = Vq + vAv, 
ot 

div v = 0. 

Noether's theorem can then be applied. The infinitesimal generators of the couple of equations 
(Navier-Stokes equations and their adjoint equations) are: 

d 

X °-di> 
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d d d d d d 

X ij = X io htt i» V i^~ ' * = l> 2 i J > h 

OXi OXj OUi OUj OVi OVj 

^ = ^t)^ + a[(t)^-x t aUt)^ p , , = 1,2,3, 
d d d d d 

Y l = 2t ^l + X ka U ka 2 P* < 7^ - > 

ot oxk out op oq 

, 9 d d 

= [XjUi XiUj)-g^ + Xj xi , i = 1, 2, j ' > i, 

X[= ( Xt a'(t)- Ul a(t))—-a(t) — , i = 1,2,3, 

where C, the c^'s, and u are arbitrary scalar functions. 

Conservation laws for Navier-Stokes equations can be deduced from these infinitesimal gene- 
rators. However, that requires non-trivial calculations and is not done in this paper. 

In the last section, we will prove that a model which is consistent with the second law of 
thermodynamics, i.e. such that the total dissipation remains positive, is stable. 

B Stability of thermodynamically consistent models 

After an eventual change of variables such that u vanishes along the boundary V of the domain 0, 
the filtered equations can be written in the following form: 

9u , „ 1 — 

— + div(w <g) u) + -Vp - div T - T s = F, 

ot p 

divTZ = (29) 

associated to the conditions 

u = sur r, 

11(0, x) = 7(2;) on 



/ p(t,x) dx = Vie [0,t f ]. 
Jn 



F is an appropriate function of t and x and tj is the final observation time. 

Proposition. Let (u,p) be a regular solution of (j29|) where T s is symmetric and verifies the 
condition: 

tr[(T - T S )S] > 0. 

Then: 



_ ft} 

\u(t,x)\\ L2{n) < \h(x)\\ L 2 (n) + / \\F(T,x)\\ L 2 (n) dr. 

Jo 



This proposition ensures a finite energy when the model conforms to the second law of 
thermodynamics. 



Consequences of Symmetries on the Analysis and Construction of Turbulence Models 



19 



Proof. Let (., .) denote the scalar product of L 2 (Q) and (u,p) a regular solution of (|29j) . From 
the first equation of (|29|) and the boundary condition, we have: 

/ du _\ _ _ 1 __ _ _ 

I — ,u 1 + b(u,u,u) - - (p,divw) + (T - T s ,Vu) = (F,u), 

where b is defined by the trilinear form 

b{u l , u 2 , u 3 ) = ( div(w x ® u 2 ), u 3 ) . 

From integrals by parts, the boundary condition and the divergence free condition, it can be 
shown that b(u, u,u) = and (p, divU) = 0. Since (T — T s ) is symmetric, it follows that 

f du \ _ 

f —,u\ +(T-T s ,S) = (F,u). 
Consequently 

Il(S,u) + (T-T.,S) = (F,S). 



Now, using the main hypothesis, we have: 
(T-T S ,S)= f tr[(T - T S )S] dx > 



and then 
1 d 



(u.u) < (F,u). 

2 dt v ; ~ v ; 



Consequently, 

d 



l w IU 2 (n)-r^ll w llL 2 (n) < (-f 1 ^) < H- F1 |lL 2 (n)ll' u llL 2 (f7)- 



After simplifying by the L 2 -norm of u and integrating over the time, it follows: 

\\u(t,x)\\ L 2 {n) < \\~/{x)\\ L 2 (n) + [ ||F(r,aj)||La(n) dr 

J o 

< ll7(«)llL 2 (n) + / \\F(T,x)\\ L 2 m dr. 



This ends the proof of the proposition. 
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